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We present measurements of the thermal conductivity n and the thermal expansion a of NdMnC>3 
and TbMnOs. In both compounds a splitting of the 4/ multiplet of the 7? 3+ ion causes Schottky 
contributions to a. In TbMnOs this contribution arises from a crystal-field splitting, while in 
NdMn03 it is due to the Nd-Mn exchange coupling. Another consequence of this coupling is a 
strongly enhanced canting of the Mn moments. The thermal conductivity is greatly suppressed in 
both compounds. The main scattering process at low temperatures is resonant scattering of phonons 
between different energy levels of the 4/ multiplets, whereas the complex 3d magnetism of the Mn 
ions is of minor importance. 

PACS numbers: 74.72.-h, 66.70.+f 



I. INTRODUCTION 

The search for magnetoelectric materials with the pos- 
sibility to influence magnetic (electric) ordering by an 
electric (magnetic) field has greatly increased the inter- 
est in so-called multiferroic materials, in which magnetic 
and ferroelectric ordering phenomena coexist 1 . The or- 
thorhombic rare-earth manganites i?MnC>3 are particu- 
larly important here, since for R = Gd, Tb, and Dy a fer- 
roelectric phase develops within a magnetically ordered 
phase. These compounds show complex magnetic struc- 
tures driven by frustration effects, and there is evidence 
that the ferroelectric order is related to a cycloidal mag- 
netic ordering. Since the different magnetic and elec- 
tric phase transitions strongly couple to lattice degrees 
of freedomS^i^, one may expect a strong influence of 
these ordering phenomena and the related low-lying ex- 
citations on the phonon thermal conductivity. In fact, 
recent zero-field thermal conductivity measurements of 
various i?Mn03 compounds by Zhou et al£ seem to sup- 
port such a view. The thermal conductivity is found 
to be drastically suppressed for R= Tb and Dy, while 
it shows a rather conventional behavior for the other 
i?MnC>3 compounds. Thus, the authors of Ref. H inter- 
preted the suppressed thermal conductivity of R= Tb 
and Dy as a consequence of the complex ordering phe- 
nomena in these compounds. Based on the above con- 
jecture we have studied the magnetic-field influence on 
the thermal conductivity of multiferroic TbMnC>3 and 
of NdMnC>3 which shows a more conventional antiferro- 
magnetic order. A detailed analysis of the thermal con- 
ductivity in combination with results from thermal ex- 



pansion measurements of NdMnC>3 and TbMnC>3 reveals 
that for both compounds the suppression of the thermal 
conductivity is largely determined by resonant scattering 
of phonons by different energy levels of the 4/ orbitals. 
Therefore, our new results are in strong contrast to the 
interpretation proposed in Ref. [^. 

The starting compound of the i?MnC>3 series, 
LaMn03, crystallizes in an orthorhombic crystal struc- 
ture (Pbnm) with a GdFeC>3 type distortion. If La is 
replaced by smaller rare-earth ions the GdFeC>3 type dis- 
tortion increases, which causes a decreasing Mn-O-Mn 
bond angle. In LaMn03, a Jahn- Teller ordered stated is 
realized below Xjt ~ 750 K, and an A-type antiferromag- 
netic (AFM) ordering of the Mn moments develops below 
T^f n m 140 K. This type of ordering is characterized by a 
ferromagnetic alignment of the magnetic moments within 
the ab planes, and an antiferromagnetic one along the c 
axisZ&2. A Dzyaloshinski-Moriya (DM) type interaction 
Jdm causes a canting of the spins towards the c direction 
resulting in a weak ferromagnetic moment (Mwf)- This 
A-type AFM ordering remains for R = Pr,. . . ,Eu, but the 
Neel temperature is successively suppressed. There are 
three main exchange couplings between the Mn moments: 
the ferromagnetic nearest-neighbor (NN) coupling (J™) 
within the ab plane, the next-nearest neighbor (NNN) 
antiferromagnetic exchange interaction («7i| M ) within 
the ab plane, and the antiferromagnetic NN interaction 
J^ FM along the c direction. A larger distortion, i.e. a de- 
creasing Mn-O-Mn bond angle, suppresses J^ M , whereas 
j AFM hardly change o 10 d 1 d 2 . The increasing frustration 
between J FM and J,f FM destabilizes the A-type AFM 
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ordering and finally leads to complex ordering phenom- 
ena for R= Gd, Tb, and Dyl^WVVVVVQ, For 
R = Dy. . .Lu, i?Mn03 crystallizes either in a GdFeOv 
type or in a hexagonal structure depending on the growth 
technique, while for R = Er, . . . , Lu usually the hexago- 
nal structure is realized 21,22 . 

The presentation of our results in the subsequent sec- 
tions is organized as follows. After a description of the 
experimental setup we first concentrate on the zero-field 
data obtained on NdMnC>3. Then we discuss the influ- 
ence of a magnetic field applied cither along the c axis or 
within the ab plane of NdMnC>3. In the last subsection 
our results obtained on TbMnC>3 will be analyzed. 



II. RESULTS AND DISCUSSION 



A. Experimental 



The NdMnC>3 single crystal used in this study is a 
cuboid of dimensions 1.65 x 1.85 x 1.2 mm 3 along the a, 
b, and c direction, respectively. It was cut from a larger 
crystal grown by floating-zone melting 2 ^. Magnetization 
and specific heat data of the same crystal are reported 
in Ref. [24. The measurements on TbMnC>3 have been 
performed on different small single crystals which were 
cut from a larger crystal grown by floating-zone melt- 
ing in an image furnace 14 . The thermal conductivity 
was measured by a standard steady-state technique us- 
ing a differential Chromel-Au+0.07%Fe-thermocouple 25 . 
For NdMnC>3, we studied Kb with a heat current along 
the b axis in magnetic fields applied either along a, 6, 
or c. All these measurements were performed with one 
set of heat contacts using either a 140 kOe longitudinal- 
field or a 80 kOe transverse-field cryostat. For TbMn03, 
we present measurements of K a , i-e. with a heat current 
along the a axis. Here, we used different configurations 
in order to allow measurements up to 140 kOe for all 
three field directions. Unfortunately, the TbMnC>3 crys- 
tal cracked when we increased the field above 110 kOe for 
H || c. This problem also occurred on another TbMnC>3 
crystal during our thermal-expansion measurements for 
the same field direction and we suspect that the crys- 
tals break because of strong internal torque effects^. The 
longitudinal thermal expansion coefficients on have been 
measured along all three crystal axes i = a, b, and 
c using different home-built high-resolution capacitance 
dilatometer a 26 ' 27 ! 28 . For the field-dependent measure- 
ments, we concentrate on ab of NdMnOs with H \\ c and 
H || b up to maximum fields of 140 kOe. Measurements 
with H || a were not possible due to large torque effects. 
For TbMnC>3, we only present zero-field data of cti, since 
the field influence is discussed in detail in Ref. \4. 



B. NdMnO :i in Zero Magnetic Field 

Figure QJi shows the uniaxial thermal expansion coeffi- 
cients cti (i = a, b, c) of NdMnC>3. The Neel transition at 
T^f n ~ 85 K causes large anomalies along all three crys- 
tallographic axes. The sign of the anomaly is positive for 
ab and a c , while it is negative for a a . From Ehrenfest's 
relations it follows that the sign of the anomaly of ai cor- 
responds to the sign of the uniaxial pressure dependence 
of Tn. This means, for example, T^f n would increase un- 
der uniaxial pressure applied either along the b or c axis. 
The different signs of the uniaxial pressure dependencies 
of Tn as well the suppression of Tn in the i?Mn03 series 
with increasing Mn-O-Mn bond angle can be essentially 
traced back to the frustration between the NN J™ and 

the NNN Jjj^ FM in the ab planes. For a more detailed 
discussion we refer to Rcfs. EH 

Figure [TJ) shows the zero-field thermal conductivity Kb 
(left scale) together with the specific heat (right scale; 
data from Ref. |24| ). There is a A-like peak in the spe- 
cific heat at T^f n ~ 85 K and at the same temperature Kb 
has a sharp minimum. Above Tn , Kb increases monotoni- 
cally, which contradicts conventional phononic behavior. 
Below Tn, k has a maximum at 25 K with a rather low 
value of 7W/Km. Around 8K the specific heat shows 
a Schottky peak, which arises from a splitting of the 4/ 
ground-state doublet of the Nd 3+ ions^.. In general, a 
two-level Schottky peak is described by 
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r 1 r 2 exp(-A/T) 
(r 1+ r 2 exp(-A/T)) 2 
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where tj and r 2 are the degeneracies of the involved lev- 
els and A the energy splitting. For NdMnC>3 we obtain 
A ~ 21 K and t\ = t 2 = 1. The 4 ig/2 ground-state 
multiplet of a free Nd 3+ ion is ten-fold degenerate, but 
splits to five doublets in the orthorhombic crystal field 
(CF). To our knowledge, no neutron scattering inves- 
tigations of the CF splitting of NdMnOs are available. 
However, in the related compounds Nd^403 with A =Ni, 
Fe, and Ga a splitting of the order of 200 K between the 
ground-state doublet and the first excited doublet has 
been measure d 29 ' 30 i 31 ' 32 . Since a similar CF splitting is 
expected for NdMnOs, the observed Schottky anomaly 
cannot arise from a thermal population of the first excited 
doublet. Instead it has to be attributed to a zero-field 
splitting of the Nd 3+ ground-state doublet, which arises 
from the exchange interaction between the canted Mn 
moments and the Nd moments^. A Schottky peak can 
also occur in ai and its magnitude and sign are given by 
the uniaxial pressure dependence of the energy garj 4 * 3 ^. 
In NdMnOs, an obvious Schottky contribution is only 
present for ab, which will be analyzed in detail below. 

In an insulator the heat is transported by phonons. 
This can be described by ft oc Gv s t, where C is the spe- 
cific heat, v s the sound velocity, and I the mean free path 
of the phonons. At low temperatures I is determined by 
boundary scattering, and k follows the T 3 dependence of 
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FIG. 1: (Color online) a) Thermal expansion of NdMnOa 
along all three crystallographic axes, b) Thermal conductivity 
measured with a heat current along the b axis (left scale) and 
specific heat (right scale; data from Ref. I24T) of NdMnOa. 



the specific heat. At intermediate temperatures k tra- 
verses a maximum with a height strongly determined by 
scattering of phonons by defects. At high temperatures 
C approaches a constant, and k roughly follows a 1/T 
dependence due to Umklapp scattering. 

Figure^ compares the zero-field thermal conductivity 
of NdMn0 3 and NdGaO^. In NdGa0 3 , the expected 
temperature dependence of conventional phononic ther- 
mal conductivity is observed, and in the whole temper- 
ature range K is significantly larger than k of NdMn03. 
This difference shows that additional scattering mecha- 
nisms are acting in NdMnOa. Fot a quantitative analysis 
we use an extended Debye model 3 -, which yields 
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e D /T 



t(x,T) 



(e* - l) z 



T dx. 



(2) 



Here, 8d is the Debye temperature, v s the sound veloc- 
ity, to the phonon frequency, x — huj/ksT, and r(x, T) 
the phonon relaxation time. The scattering rates of dif- 
ferent scattering mechanisms, which are independent of 
each other, sum up to a total scattering rate 
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The four terms on the right-hand side refer to the 
scattering rates on boundaries, on planar defects, on 
point defects, and to phonon-phonon Umklapp scatter- 
ing, respectively. The mean free path cannot become 
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FIG. 2: (Color online) a) Zero-field thermal conductivity of 
NdMnOs (O) and NdGa03 (A) on a logarithmic « scale. The 
lines are calculated within the Debye model using the pa- 
rameters D = 600 K, v s = 6000 m/s, P = 3.9 ■ I0" 43 s 3 , 
U = 6 ■ I0~ 31 s/K, u = 5.6, and l min = 9 A (see text). For 
the solid and dashed lines magnetic scattering has been in- 
cluded (e = 4.5 • I0" 41 m 3 s 3 /KJ) and excluded (e = 0), re- 
spectively, b) Specific heat of NdMn0 3 (o, Ref. H3) and of 
NdGa03 (A, Ref. I34T ) as a reference compound. The mag- 
netic contribution C m (T) (?) due to Mn spin excitations 
is estimated by the difference of both curves. Inset: Mag- 
netic entropy S m = J C m (T) /T dT (solid line) and expected 
Sm = N A k B ln(2S + l) ~ 13.4 J/molK (dashed) for the S = 2 
spins ofMn 3+ . 



smaller than the lattice spacing giving a lower limit for 
k. This is taken into account in equation ^ by a 
minimum mean free path £ m j n and replacing t(x, T) by 
max{rs(x, T), £ m i n /v s } (Ref. 36). Theoretical investiga- 
tions of the scattering of phonons by magnetic excitations 
yield an additional scattering rat o 37 ' 38 



eT 2 C m (7> 4 , 



(4) 



where e describes the scattering strength, and C m is the 
magnetic contribution to the specific heat. Note, that 
fluctuations may cause a sizeable C m also above Tn- 

In order to determine C m the other contributions to C 
have to be subtracted. In NdMnOa this background con- 
tribution Cb g arises from acoustic and optical phonons as 
well as the Schottky specific heat of the 4/ CF excitations 
of Nd 3+ . Since a calculation of Cb g with the required pre- 
cision is not possible, we use NdGaOa as a reference com- 
pound. Due to the structural similarit y 39 ' 40 , the phonon 
spectrum and the CF splitting are presumably very sim- 
ilar in NdGaOs and NdMnOa, apart from the additional 
splitting of the CF doublets of NdMnOa ■ Figuredb shows 
the specific heat of NdMnO^ and NdGaO^. At high 
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temperatures the curves are indeed nearly identical. We 
estimate 

C m (T) = CNdMnOa (I") _ CNdGa0 3 (T) - C^(T). (5) 

The last term describes the Schottky contribution due to 
the splitting of the ground-state doublet and is calculated 
by equation ([1]) with A = 21 K and t\ = r% = 1. Note 
that the splittings of the excited doublets do not need 
to be considered, since their population sets in at higher 
temperature, and and as long as the splittings are not 
too large they hardly change the specific heat. In Fig. [2)3 
we display the resulting C rn , which exhibits the A peak 
at Tn and then slowly decays for T > Tn. As a test of 
our analysis, we also calculate the magnetic entropy 



Sm(T) = Io^ dT ' 



(6) 



As shown in the Inset of Fig. [^b, the experimental 
S m (T) approaches the expected ATa^e ln(25 + 1) — 
13.4 J/moleK of the Mn moments, where Aa and fcs de- 
note Avogadro's number and Boltzmann's constant, re- 
spectively. 

The dashed line of Fig. [2jt shows the thermal con- 
ductivity calculated for NdGaC>3 with the parameters 
given in the figure caption. The Debye temperature and 
the sound velocity have been estimated from the mea- 
sured specific healj^l and the other parameters have been 
adapted to fit the data. The calculated curve reproduces 
the general behavior of K well. In order to describe the 
thermal conductivity of NdMnC>3 the additional mag- 
netic scattering rate r" 1 is switched on by adjusting the 
parameter e, while keeping all the other parameters fixed. 
This calculation (solid line) describes the temperature 
dependence of k above Tn very well. However, the cal- 
culation overestimates the minimum at Tn, and it shows 
a pronounced low-temperature maximum which is not 
present in the data. In principle, the latter difference 
could arise entirely from different point defect scattering 
in NdMnC>3 and NdGaC>3. However, our magnetic-field 
dependent measurements will show that this difference 
arises to a large extent from an additional phonon scat- 
tering on the CF levels. 



C. NdMn0 3 in a Magnetic Field H \\ c 

The low-temperature Schottky contribution to the 
thermal expansion of NdMnC>3 is most pronounced for 
ab] see Fig. Q] For a two-level system this contribution 
follows from a Griineisen scaling between a and the spe- 
cific heat 3 - 3 -, and for two singlets it gives 



k B dln(A) / A 
d Pi It 



-A/T 



(l + e- A / T ) 2 
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Here, V uc is volume per formula unit. The magnitude 
and the shape of cn S ch,i(T) are entirely determined by 
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FIG. 3: (Color online) a) Thermal expansion otb (symbols) 
of NdMnC>3 for various magnetic fields H || c together with 
Schottky fits (solid lines). The arrow indicates the evolution 
of the Schottky peak with increasing field. Note the logarith- 
mic temperature scale. Inset: Energy splitting A (symbols) of 
the 4/ ground-state doublet of Nd 3 as a function of magnetic 
field. The line is a linear fit of A(H). b) Thermal conductiv- 
ity of NdMn03 for various H \\ c as a function of temperature 
on double-logarithmic scales. Inset: Sketch of the canting of 
the Mn spins along c (see text). 



the energy gap A and its uniaxial pressure dependence 



ain(A) 



. We conclude that 



gln(A) 



is rather small for uni- 



axial pressure applied along a or c, since the Schottky 
contributions of a a and a c are so small that they are 
almost entirely masked by the respective phononic con- 
tributions. In contrast, ah is clearly dominated by the 
Schottky contribution a S ch.6 up to f=a 20 K. Thus, ctb al- 
lows for a detailed analysis of the magnetic-field depen- 
dence of the splitting of the Nd 3+ ground-state doublet. 
Figure [3k. shows a b in magnetic fields up to 140 kOe ap- 
plied along the c direction. With increasing field, the 
Schottky peak shifts monotonically to higher tempera- 
tures. The solid lines in Fig.[3h. are fits via equation Q, 
which very well reproduce the experimental data and al- 
low to derive the energy gap as a function of magnetic 
field. As shown in Fig. [3)d we find a linear increase 
A(H) = 21 K + 2.25 • 10~ 4 K/Oe • H. The zero-field 
value nicely agrees with that obtained from the specific 
heat^i. The field dependence of A can be understood as 
follows. In zero field, Ao solely arises from the Nd-Mn 
exchange, which is proportional to the zero-field M^ F - 
A field applied along c increases A, on the one hand, due 
to the additional Zeeman splitting of the Nd 3+ ground- 
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state doublet. One the other hand, the canting of the 
Mn moments also increases with field, yielding an addi- 
tional increase of the Nd-Mn exchange. Thus, for H || c 
we obtain 

A(H) = ±(M^ + XMn H) + . (8) 

Here, 5 is the proportionality constant between Mwf and 
the Nd-Mn exchange, XMn the field-dependence of M WF , 
and the g factor of the Nd 3+ ground-state doublet. 

For LaMn0 3 and PrMn0 3 , values of M^ F ~ 0.1^ B 
due to the DM interaction are reporte d 24 ' . A similar 
DM interaction can be expected for NdMn0 3 , but due to 
the Nd-Mn exchange interaction additional energy can be 
gained from an enhanced splitting of the Nd 3+ ground- 
state doublet by increasing the canting of the Mn spins. 
In order to estimate this effect we calculate the single-site 
energy of a Mn 3+ ion as a function of the canting angle 5; 
see Fig.[3ji. In a first step, we consider the zero-field case 
of PrMn0 3 , where only J^ FM , Jum, and the single- ion 
anisotropy D of Mn 3+ have to be considered. Based on 
the Hamiltonian of Ref. [TJ we obtain 

E^(5) = -4jf FM S ,2 cos(2<5) - 2J DM sin5 - D[S cos S} 2 

(9) 

where S — 2 is the Mn spin. Using J^ FM = 7 K and 
D = 0.9 K from Ref. |H and S = 1.4° from M^ F = 
4/XBsin<5o = 0.1/xb^, we obtain Jdm = 5.7 K from the 
minimization condition dE/dS = at So = 1.4°. In the 
next step we include the additional energy gain in the 
Nd 3+ ground-state doublet, which is given by A(H)/2 
from equation (J5J) , and the potential energy of the Mn 
moment in a finite field H \\ c. It is reasonable to assume 
that J^ FM , Jdm, and D do not change from PrMn0 3 to 
NdMn0 3 (see also Refs. HHIIHil). Thus, we keep these 
parameters fixed and obtain 

E Nd (S,H) = E**(S) ^ - MwfWJ? . (10) 

2 kb 

The determination of the remaining parameters is 
straightforward. First dE(S,H = Q)/dS = is solved 
for H = under the additional condition that the 
zero-field value Ao = 21 K is reproduced. This yields 
M°, F = 0.65 /ib, — 9.4°, and h = 470 kOe. Then xm„ 
and (?Nd follow from a minimization of equation (|10p for 
finite fields under the condition that the observed field de- 
pendence dA/dH = 2.25 • 10~ 4 K/Oe is fulfilled. The re- 
sulting values are <7Nd = 2.2 and XMn = 2.4 • 10~ 6 ^/Oe. 
Remarkably, the weak fm moment Af^, F — 0.65 ^b of 
NdMn0 3 is strongly enhanced compared to the values 
of ~ 0.1 /is of PrMn0 3 or LaMn0 3 . This enhancement 
should be clearly visible in a magnetic structure deter- 
mination, which would be a good test of our analysis. 

Figure [3Jd shows the thermal conductivity of NdMn0 3 
up to H = 140 kOe applied along the c direction. At 
5 K the thermal conductivity increases almost linearly 
with field up to Kb — 20 W/Km. This strong field depen- 
dence weakens with increasing temperature, and around 



25 K the field dependence even changes sign and remains 
negative up T > Jn- Our data suggest that the low- 
temperature behavior of Kb arises from resonant scatter- 
ing of phonons between different levels of the 4/ multi- 
plet of Nd 3+ which causes a suppression of the phonon 
heat transport in a certain temperature range. Such a 
suppression of k by resonant scattering on 4/ states is 
well-known from the literature^ 3 -. The idea is that a 
phonon with an energy equal to the energy splitting of 
two 4/ levels is first absorbed and then reemitted. Since 
the momenta of the incoming and reemitted phonons 
have arbitrary directions, an additional heat resistance 
is caused (for more details see e.g. Ref. H3). The com- 
parison of the zero-field thermal expansion and the zero- 
field thermal conductivity data gives clear evidence that 
resonant scattering between the two levels of the split 
ground-state doublet is the cause for the strong suppres- 
sion of k at low temperatures. Since the splitting of the 
ground-state doublet increases with increasing field the 
scattering probability of the low-energy phonons system- 
atically decreases resulting in a strong increase of the 
low-temperature thermal conductivity. In the tempera- 
ture range above ~ 25 K, the situation is more complex 
as will be discussed in the following subsection. 



D. NdMnO ;i in Magnetic Fields H || a and \\b 

Figure [4^i shows the temperature-dependent thermal 
expansion ab for magnetic fields H || b up to 80kOe. 
Here, the behavior of the Schottky contribution is dif- 
ferent to that observed for H \\ c. For small fields (H < 
20kOe) almost no effect occurs while for higher fields the 
peak height continuously decreases until it disappears 
completely for H ~ 80 kOe. The maximum of the peak 
weakly shifts to higher temperature when the field is in- 
creased from to 60kOe. This weak increase is a con- 
sequence of the perpendicular orientation of the external 
magnetic field with respect to the exchange field arising 
from M\vf||c Thus, the total effective field is given by 
the vector sum of both contributions, which for small ex- 
ternal fields only weakly increases. The decrease of the 
peak height and its disappearance at 80 kOe suggest that 
the pressure dependence of the energy gap also decreases 
with field and finally vanishes. Whether this is really 
the case is, however, not clear because around 100 kOe a 
spin- flop transition takes place for this field direction 24 . 
Thus, different energy scales have to be considered, which 
prevent a simple analysis of the pressure dependencies 
via a Griineisen scaling 33 . As displayed in Fig. 0b, ab 
for H > 100 kOe| 1 6 has another anomaly, which we at- 
tribute to the spin-flop transition. This anomaly strongly 
shifts towards higher temperatures with further increas- 
ing field, i.e. the phase with M\vf||c becomes less stable 
towards lower temperatures. 

Figure [5] shows k for H \\ a and which is suppressed 
in the entire range T < for both field directions. 
This behavior is in clear contrast to the strong low- 
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temperature increase of n for H || c, whereas the weak 
decrease of k above about 30 K is rather similar for all 
three field directions. Thus we conclude that the field 
dependence of Kb is determined by two different scatter- 
ing mechanisms. Firstly, there is the resonant scattering 
within the split ground-state doublet. This scattering 
is strongly suppressed at low temperature for H || c be- 
cause of the increasing splitting. As discussed above, the 
splitting increases much less for H \\ a and H || 6, since 
the effective field increases only weakly for H ± Mwf- 
Nevertheless, a low-temperature increase of Kb should 
occur for H \\ a and H || b if this resonant scattering was 
the only field dependent scattering process. To explain 
the observed decrease of Kb requires another scattering 
process which increases with magnetic field for all three 
field directions. This second process is only visible when 
the field dependence of the first resonant one is weak, i.e., 
for H || a and H \\ b in comparatively weak fields and for 
H || c at T > 25 K. Since the field dependence of K a es- 
sentially vanishes slightly above Tn we suspect that this 
second field dependent scattering process is related to 
scattering of phonons by magnons, but the presence of 
higher-lying CF levels of Nd 3+ could also play a role. 

E. TbMn0 3 

TbMn03 is the first compound of the i?MnC>3 series 
in which ferroelectricity has been established over a large 
temperature and magnetic-field range. The phase dia- 
gram was first explored by polarization and magnetiza- 



tion data in Ref. [TH and recently it has been refined by 
thermal expansion measurements^. In zero field, the sys- 
tem transforms from a paramagnetic to an incommensu- 
rate antiferromagnetic phase (HTI) at Tn — 41 K. At 
Tfe — 27 K, a transition occurs into another incom- 
mensurate antiferromagnetic phase (LTI) with a differ- 
ent propagation vector. This phase is ferroelectric with 
a polarization along c. The phase boundaries at Tn and 
Tfe hardly depend on a magnetic field, but for 7T || a 
and || b a transition from the LTI phase to a commen- 
surate antiferromagnetic phase (LTC) occurs below Tfe, 
which is accompanied by a polarization flop from P \\ c to 
P || a. The main difference between H \\ a and H \\ b is 
the much larger hysteresis of the LTI-to-LTC transition 
for H || a. A magnetic field along c causes a transition 
into a paraelectric canted AFM phase above ~ 7 T. 

Figure [S{a-c) shows the thermal conductivity of 
TbMnC>3 along the a direction. In zero field, K a has 
a broad minimum around 80 K and a weak maximum 
at T ~ 34 K with a relatively low absolute value k « 
3W/Km. Figure [BJi displays the thermal expansion of 
TbMnC>3 in zero field along the a, b, and c axes. The 
transitions at Tn and Tfe cause anomalies of a along 
all three crystallographic directions, which are discussed 
in detail in Ref. [j. In addition, for all three direc- 
tions broad Schottky contributions of different signs are 
present around 80 K, which originate from the 4/ states 
of Tb 3+ . The 7 T 6 state of the free Tb 3+ ion splits into 
13 singlets in an orthorhombic CF. To our knowledge, 
the 4/ energy level scheme of Tb 3+ in TbMn03 is not 
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T(K) 

FIG. 6: (Color online) a-c): Thermal conductivity n a of 
TbMnOa in magnetic fields along the different crystal axes, 
d) Thermal expansion at of TbMnOa with i = a, b, and c. 
The vertical line is a guide for the eyes. 



known— and we are not aware of any investigations of 
the crystal-field splitting of related TbAOs compounds. 
Thus, a detailed analysis of this Schottky contribution 
to a is not possible. As a rough estimate of the relevant 
energy scale a fit of the thermal expansion data by equa- 
tion ([7]) yields an effective energy gap of m 190 K. The 
comparison of the thermal conductivity with the thermal 
expansion data shows that the extrema of ai occur close 
to the minimum of k. This correlation strongly suggests 
that the minimum of n a at ps 90 K is caused by reso- 
nant scattering of phonons between different CF levels of 
Tb 3+ . Such an interpretation is further supported by the 
magnetic- field dependencies of K a - For H \\ b and H || c, 
K a is suppressed over a broad temperature range, whereas 
it is enhanced for H \\ a. These large field dependencies 
up to high temperatures are rather unusual, and clearly 
not related to the low-tcmperature ordering phenomena 
below 41 K. The influence of the various transitions be- 
low 41 K on K a is weak. We only observe small dips at 
the transition temperatures Tn and Tfe as indicated by 
the arrows in Fig. [5] The shape of the thermal conductiv- 
ity curves remains essentially unchanged. We conclude 



that these ordering transitions play little role for K a , since 
the low K a is completely dominated by the scattering of 
phonons by the 4/ CF levels of Tb 3+ . 

Below Tfe additional magnetic-field dependencies are 
observed. One can distinguish the measurements with 
if || a or b, where the LTI-to-LTC transition occurs, from 
those with H || c, where the system turns into the para- 
electric phase for sufficiently large fields. For the latter 
case, a sharp increase of K a is observed when the para- 
electric phase is reached, see Fig. [SJ;. The pronounced 
hysteresis of the 80 kOe curve reflects the first-order na- 
ture of this transition. Apart from the different transition 
temperatures, the curves for 80 kOe and 110 kOe are al- 
most identical. For H = 80kOe || b, n a increases at the 
LTI-to-LTC transition (w 18 K; see Fig. |6b). Although 
this field dependence is of opposite sign compared to the 
high-temperature field dependence for H || b, it is diffi- 
cult to separate both effects, because the scattering by 
the CF levels of Tb 3+ probably causes a magnetic-field 
dependence of k down to the lowest temperature. Since 
the LTI-to-LTC transition is accompanied by a polariza- 
tion flop from P || a to P \\ c, one may suspect that 
the increase of n a is related to the formation of ferro- 
electric domains. We have, however, ruled out this pos- 
sibility by measurements of the electrical polarization as 
well as of the thermal conductivity under application of 
large electrical fields (not shown). Although the domain 
formation could be clearly seen in the polarization mea- 
surements, no electric-field influence was detectable in 
n. Thus, another explanation for the suppressed thermal 
conductivity in the LTI phase is needed. Probably, it is 
the incommensurability itself, which causes an additional 
thermal resistance, because the crystal symmetry is low- 
ered in the LTI phase. Accordingly, one should expect 
a similar behavior for H || a, but the broadening of the 
LTI-to-LTC transition due to its large hysteresis^ and the 
superimposed high-temperature field dependence prevent 
a separation of the different effects at low temperatures. 
The incommensurability is also consistent with the jump 
of k for H || c, since the paraelectric phase is commensu- 
rate and therefore of higher symmetry. 



III. CONCLUSIONS 

We have studied the thermal expansion and thermal 
conductivity of NdMn03 and TbMnOa under applica- 
tion of large magnetic fields. The thermal conductivity 
of NdMnOa is a very unusual. The Neel transition at 
Tn— 85 K leads to a strong suppression of the phonon 
thermal conductivity over a large temperature range. In- 
cluding a magnetic scattering rate proportional to the 
magnetic specific heat allows us to describe the thermal 
conductivity from Tn to room temperature. At low tem- 
peratures the thermal conductivity is further suppressed 
by another scattering mechanism. The 4/ ground-state 
doublet of Nd 3+ is split (Ao ~ 21 K) by an exchange 
interaction with the canted Mn moments. Our analy- 
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sis suggests a significant enhancement of the Mn canting 
angle in NdMn03 compared to that in PrMnC>3 as a con- 
sequence of this Nd-Mn interaction. The splitting of the 
ground-state doublet thereby allows for resonant scatter- 
ing of phonons which causes the additional suppression 
of k in zero field. The analysis of the thermal expan- 
sion in magnetic fields up to 140 kOe reveals that A(-ff) 
strongly increases in magnetic fields H \\ c. This increase 
of A (H) shifts the effectiveness of the resonant scattering 
processes towards higher temperature and causes a dras- 
tic increase of n at low temperatures. For H | c, with 
increasing temperature a gradual change occurs leading 
to a suppression of n. A similar suppression is present 
for H || a and H \\ b in the entire low-temperature range. 
This requires the presence of a second field-dependent 
scattering mechanism, which may be related to scatter- 
ing of phonons either by magnons or by higher-lying CF 
levels. 

TbMn03 also exhibits a strongly suppressed thermal 
conductivity over the entire temperature range. The 
clear correlation of the temperature dependencies of k 
and of the uniaxial thermal expansion coefficients a en- 
ables us to conclude that the dominant mechanism sup- 
pressing k is resonant scattering of phonons by the 4/ CF 
levels of Tb 3+ . The interpretation of Ref. [H that the low 



absolute values of the thermal conductivity of TbMnC>3 
should be caused by the complex magnetic and electric 
ordering phenomena is ruled out by our data. In contrast, 
the complex transitions of TbMnOs only cause very weak 
anomalies in k at Tn and T FE . A somewhat larger influ- 
ence is present at the transitions induced by finite mag- 
netic fields. The LTI-to-LTC transition for H || a, b as 
well as the transition to the paraelectric phase for H \\ c 
cause an increase of the thermal conductivity. Probably, 
this increase of n arises from the incommensurability of 
the LTI phase, which is transformed to a commensurate 
phase of higher symmetry for all three field directions. 
We also found that the ferroelectric domain structure 
has no measurable influence on the heat transport in 
TbMn0 3 . 
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